Comment on "Photonic Band Gaps: Noncommuting limit and 

the 'Acoustic Band'" 



Didier Felbacq 
LASMEA UMR-CNRS 6602 

Complexe des Cezeaux 
63177 Aubiere Cedex, France 
(February 6, 2008) 

This comment concerns the homogenization of 2D dielectric photonic crystals, and the 
fact that the limits k — > (k is the Bloch vector) and e — > +oo do not commute for 
p-polarized waves (it is so for s-polarized waves, which is a straightforward case). This 
result has been claimed to be true by Nicorovici & al. in a series of paper |l|-|3j and it 
has been claimed to be false by Krokhin & al. in a comment The point of this note 
is to make the situation clear one for all, that is to give a mathematically clean derivation 
of the result by Nicorovici and to prove that it is right. The first point is to use a clear 
definition of a homogenization process: rather than letting tend the Bloch vector to zero it 
is in my opinion clearer to deal with a finite-size photonic crystal, contained in a bounded 
domain Q, with period r\ (the period is a contracted cell r]Y, where Y = [0,1[ 2 and 8 is 
the filling ratio in Y see fig. 1 for notations) and a fixed wavenumber ko, in which case 
for an incident field u l the total field satisfies in p-polarization div(e~ 1 Vtt r? ) + k^u^ = 0, 
where e v represents the relative permittivity of the rods which are homogeneous circular 
cross-section rods, the permittivity of one rod being equal to e s . Then we study the limit 
of u n when rj — > (in case of an infinite crystal with no incident field, this amounts to let 
k tend to zero for a Bloch wave). I have shown in a previous paper H that u v tends to uq 
satisfying div(ef M 1 m Vwo) + k^u = where 
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and (f) e is a term defined in M. If we let formally tend e s to infinity we get 
v-i ± fdw\ 

K dyi/ 

mean, of the following problem 



1 — 9 + 0oo) where <^oo = ( — — ) , w being the unique y-periodic solution, with null 
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and {.) Y denotes averaging over Y. 

Now let us deal directly with the infinitely conducting crystal. At step rj the field satisfies 

Ou 

Au v + k^Urf = in the complementary of the rods, which are denoted by T v , and = 
on dTjj. Our result is 
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Theorem 

When rj tends to 0, u v tends to u satisfying div(e^ o 1 VM ) + k 2 u = where 

= f (1 - 9) (1 - 9 + infi 
\ 1 outside 

A simple comparison with ([I]) leads to the evident conclusion that 
Corollary 

The limits n — > and £ — > +oo do not commute. 
Proof 

The field it,, is null inside the rods but we can define a function u v such that u v = u v 
outside T v and u v is in the Sobolev space H\ oc (M 2 ), so that (l — 1^) Using now a 

test function in the Schwartz space V (ft), we find 

- / Vu v V(f)d 2 x + k 2 Q I (l - l r J u v (f)d 2 x = 
Jn Jn 



Assuming that (it,,) is bounded in L 2 (Q) it can easily be shown that (u v ) is bounded in 
H 1 (f2) because it satisfies a standard Helmholtz equation Au v + k 2 (l — 1t v ) u v — . Then, 
up to the extraction of a subsequence we have u v — > uq strongly in 1? (Q), u = (1 — 9) u 
and Xrj = Vn^ — ^ Xo weakly in H 1 ($7) so that: — J n Xo d 2 £ + fcg J n u (p d 2 x = meaning 
that 

div (xo) + k 2 u = 0. 

We then have to find an expression for xo- We set iOj = w v + Xi, where w v = w (note 

that Wi — ^ Xi), we have — J Q Vn^V (4>w i) d 2 x + k 2 f n u v (f)Wid 2 x = — J n XriV(f)Wid 2 x + 
J Q V^VuiiUr, d 2 x + k\ f n u v (/)Wi d 2 x = 0. Then letting n tend to 0, we obtain 

Jn Jn 

this shows that Xo- e « — ((Vif) y + (1 — 9) e,) Vmq. The theorem and corollary follow by the 
rotational invariance of the problem. 

As a conclusion, we might suggest that a good mathematical background should make 
it possible to avoid any polemical discussions over these issues. 
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